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Abstract 

DNA sequences were obtained from the mitochondrial cytochrome-Z? gene of nine speci¬ 
mens of Nyctomys sumichrasti collected in Mexico and Central America. Phylogenetic analysis 
(Bayesian Inference) of these sequences document heretofore unrecognized patterns in genetic 
diversity among phylogroups that: 1) indicated substantial levels of genetic divergence among 
phylogroups; 2) resulted in paraphyly of taxa currently recognized as N. sumichrasti ; and 3) 
argued for a re-assessment of the current taxonomy of Nyctomys and perhaps recognition of a 
new genus. 
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Introduction 


Nyctomys sumichrasti (Sumichrasfs vesper rat; 
De Saussure 1860) is an arboreal rodent (Cricetidae, 
Tylomyinae, Musser and Carleton 2005) distributed 
from Jalisco, Mexico, to Panama (Hall 1981; Reid 
2009). It inhabits evergreen lowlands and lower mon¬ 
tane regions including cloud, secondary, riparian, and 
semi-deciduous forests (Hall 1981; Sanchez-Hemandez 
et al. 1999; Cervantes et al. 2004; Hunt et al. 2004). 
Typically, N. sumichrasti prefers middle and upper level 
forest strata, rarely descending to the ground (Emmons 
1997; Timm and LaVal 2000; Hunt et al. 2004). Given 
its arboreal to semi-arboreal behavior, specimens of A. 
sumichrasti are relatively rare in most museum collec¬ 
tions and consequently a paucity of information is avail¬ 
able concerning its ecology, genetic and morphologic 
variation, and systematics. 


Nine subspecies of Sumichrasfs vesper rat are 
recognized ( colimensis, costaricensis, decolorus, 
florencei, nitellinus, pallidulus, salvini, sumichrasti, 
and venustulus ; Hall 1981), however the majority of 
information is restricted to the original description of 
each taxon and little has been reported concerning the 
taxonomy, geographic and genetic variation, or relation¬ 
ships among subspecies. Musser and Carleton (2005) 
noted that specimens in the United States National 
Museum could be separated into two groups based on 
variation in the carotid circulatory pattern (complete 
versus derived) and molar-root number (three versus 
four). Their groups corresponded to populations oc¬ 
curring north and west of the Isthmus of Tehuantepec 
(pallidulus and sumichrasti ) versus those to the south 
and east (< costaricensis, decolorus, florencei, nitellinus. 
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and venustulus). No comments were made concern¬ 
ing colimensis and salvini. Karyotype data have been 
reported for two subspecies, with N. s. colimensis be¬ 
ing polymorphic (2n = 50-52, FN = 52-54; Lee and 
Elder 1977; Haiduk et al. 1988) and N. s. florencei 
possessing a 2n = 50 and FN = 52 karyotype (Bradley 
and Ensink 1987). 


Herein we assess genetic variation among nine 
specimens representing five of the nine subspecies 
(Fig. 1) of TV. sumichrasti. To accomplish this, DNA 
sequences from the mitochondrial cytochrome-Z? gene 
(Cytb) were obtained and then evaluated in a phylo¬ 
genetic context. 


100° w 90° w 80° w 



Figure 1. Map depicting collecting localities for samples examined in this study. Subspecies and distributions are 
recognized following Hall (1981). Numbers indicate localities as referred to in the Appendix. 


Methods 


DNA sequences. —Mitochondrial DNA was iso¬ 
lated from approximately 0.1 g of frozen liver tissue 
using the Gentra Puregene Cell and Tissue Kit (Gentra 
Systems, Minneapolis, Minnesota) or from 0.1 g of liver 
preserved in 95% ethanol using the Qiagen DNeasy™ 
Tissue Kit (Qiagen Inc., Valencia, California). For 
most specimens, the entire Cytb gene (1,143 bp) was 
amplified by polymerase chain reaction (PCR, Saiki et 
al. 1988) using the following primers: MVZ05 (Smith 


and Patton 1993) and PER03 5 (Tiemann-Boege et al. 
2000), or with primer pairs as follows: L14724 (Irwin 
et al. 1991) with CB3H(Palumbi 1996), andFl (Whit¬ 
ing et al. 2003) with MVZ14 (Smith and Patton 1993). 
For other specimens, an approximately 400 bp fragment 
(at the 5’ end) was amplified using primers MVZ05 and 
400R (Peppers and Bradley 2000). The thermal profile 
for PCR reactions was as follows: initial denaturation 
at 95°C for 2 min, followed by 35 cycles of denatur- 
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ation at 95°C for 1 min, annealing at 51°C for 1 min, 
and extension at 72°C for 2 min, with a final extension 
at 72°C for 7 min. PCR products were purified using 
ExoSAP-IT (USB Products, Cleveland, Ohio). Prim¬ 
ers used to cycle sequence the products consisted of 
MVZ05 and MVZ14 (Smith and Patton 1993), PER03 5 
(Tiemann-Boege et al. 2000), 870R (Peppers et al. 
2002), FI (Whiting et al. 2003), HI5149, and L14841 
(Irwin et al. 1991). Cycle sequencing reactions were 
purified using isopropanol cleanup protocols. Purified 
products were sequenced with an ABI 3100-Avant 
automated sequencer or with an ABI 377 automated 
sequencer using ABI Prism Big Dye version 3.1 ter¬ 
minator technology (Applied Biosystems, Foster City, 
California). Resulting sequences were subsequently 
assembled and proofed using Sequencher 4.9 software 
(Gene Codes, Ann Arbor, Michigan); chromatograms 
were examined to verify all base changes. The resulting 
DNA sequences were deposited in GenBank and the 
accession numbers are listed in the Appendix. 

Data analyses. —Two approaches were used for 
data analysis. First, representatives of all tribes con¬ 
tained in the Neotominae (Bradley et al. 2004; Reeder 
and Bradley 2004, 2007; Musser and Carleton 2005; 
Reeder et al. 2006; Miller and Engstrom 2008) were 
included to test for monophyly of the subfamily Tylo- 
myinae (Musser and Carleton 2005). The taxonomic 
sampling for this analysis included three members of 
the Neotomini ( Hodomys alleni , Neotoma mexicana , 
and Xenomys nelsoni ), one member of the Ochro- 
tomyini ( Ochrotomys nutalli ), two members of the 
Baiomyini (Baiomys taylori and Scotinomys teguina), 
two members of the Reithrodontomyini (Isthmomys 
pirrensis and Reithrodontomys fulvescens ), and eight 
members of the Peromyscini ( Habromys lepturus , 
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Megadontomys thomasi , Neotomodon alstoni , Onycho- 
mys arenicola , Osgoodomys bander anus, Peromyscus 
californicus , Peromyscus maniculatus , and Podomys 
floridanus). For these taxa, sequences were obtained 
from GenBank (accession numbers are provided in the 
Appendix). Two members of the subfamily Sigmodon- 
tinae (Oryzomys palustris and Sigmodon hispidus) were 
used as outgroup taxa for this analysis. 

In the second analysis, and based on results from 
the first analysis, the dataset was reduced to include 
only members of the Tylomyinae ( Nyctomys , n = 9; 
Otonyctomys , n = 1; Ototylomys , n = 1; and Tylomys , 
n = 1). Given that Ototylomys and Tylomys were sister 
taxa in the first analysis, Tylomys was selected for the 
outgroup taxon. Sequences were deposited in GenBank 
(accession numbers are provided in the Appendix). 
Locality information for the specimens of Nyctomys 
is provided in the Appendix. 

The GTR+I+G model was identified by the 
Akaike information criterion in MODELTEST (Posada 
and Crandall 1998) as the most appropriate model of 
DNA evolution for both datasets. A Bayesian model 
(MrBayes; Huelsenbeck and Ronquist 2001) was used 
to obtain a phylogenetic tree and to generate support 
values (clade probabilities). The GTR+I+G model 
parameters included: a site-specific gamma distribu¬ 
tion, four Markov-chains, 10 million generations, and 
sample frequency = every 1,000th generation. After 
a visual inspection of likelihood scores, the first 1,000 
trees were discarded and a consensus tree (50% major¬ 
ity rule) was constructed from the remaining trees. The 
Kimura two-parameter model of evolution (Kimura 
1980) was used to obtain genetic distances and estimate 
divergence times between taxa. 


Results 


The first Bayesian analysis produced a tree 
(not shown) depicting a monophyletic Tylomyinae 
( Nyctomys , Otonyctomys , Ototylomys , and Tylomys) 
and Neotominae (15 genera), with each clade being 
supported with probability values = 1.00. Based on 
these results, a second Bayesian analysis (Fig. 2) was 
conducted using Tylomys as the outgroup taxon (as 
explained in the Methods) in order to better assign 


character polarity to the various sequences of Nycto¬ 
mys. This analysis produced two major clades (I and 
II) and six minor clades (A-F) that were supported 
with probability values = 1.00. Clade I was comprised 
of two groups: Otonyctomys hatti, and seven of the 
nine samples of Nyctomys. Clade II consisted of the 
remaining two samples of N. sumichrasti. 
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Tylomys nudicaudus 


Ototylomys phyllotis 


I 0-99 


1.00 


II 


B 


1.0 


Otonyctomys hatti 


c _Loo 


1.00 


D 


l.OOsi 

F 


, N. s. decolorus TTU115406 Guatemala 


^ N. s. decolorus TTU84484 Honduras 
L N. s. florencei TTU111528 El Salvador 
N. s. florencei TTU 111 529 El Salvador 
N. s. florencei TTU105094 Nicaragua 
N. s. florencei TTU100336 Nicaragua 
,N. s. colimensis TTU37743 Jalisco, MX 


, N. s. salvini BYU14466 Chiapas, MX 


N. s. sumichrasti CMC101 Oaxaca, MX 


F igure 2. Phylogenetic tree obtained from the Bayesian analysis of mitochondrial cytochrome-/? gene sequences. 
Values above branches indicate clade probability support values (only values > 0.95 are shown). 


Kimura two-parameter distance values (Kimura 
1980) were obtained for comparisons among samples 
of Nyctomys, samples of Nyctomys versus Otonyctomys, 
and between other sister genera within the Tylomyinae 
and Neotominae (Table 1). Genetic distances obtained 
from within clade comparisons of samples of Nyctomys 
ranged from 1.5% (clade F) to 10.2% (clade II), whereas 
comparisons between clades containing samples of 
Nyctomys ranged from 8.4% (clades E and F) to 20.0% 


(clades I and II). Using a sequence divergence estimate 
of 3.5% per million years for Cytb sequences (Arbogast 
and Slowinski 1998), divergence times and confidence 
intervals (obtained using the “POISSCI” function in 
MATLAB® version 6.5, The Math Works, Inc. 2006) 
within clades ranged from 0.4 ± 0.34 mya (clade F) to 
2.9 + 0.58 mya (clade II) and between clades from 2.4 
± 0.54 mya (clades E and F) to 5.7 ± 0.86 mya (clades 
I and II) (Table 1). 
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Table 1. Genetic distances estimated using the Kimura 2-parameter model of evolution (Kimura 1980) for 
selected comparisons o/Nyctomys and other members of the Tylomyinae and Neotominae. 


Comparison 

Genetic Distance 

Time Since Divergence 

Within Clades containing Nyctomys 

Clade I 

9.3% 

2.7 mya ± 0.56 mya 

Clade II 

10.2% 

2.9 mya ± 0.58 mya 

Clade B 

6.5% 

1.9 mya ± 0.49 mya 

Clade C 

3.8% 

1.1 mya ± 0.41 mya 

Clade F 

1.5% 

0.4 mya ± 0.34 mya 

Between Clades containing Nyctomys 

Clades I and II 

20.0% 

5.7 mya ± 0.86 mya 

Clades A and B 

17.7% 

5.1 mya ± 0.80 mya 

Clades C and D 

13.3% 

3.8 mya ± 0.67 mya 

Clades D and E 

14.4% 

4.1 mya ± 0.70 mya 

Clades D and F 

13.1% 

3.7 mya ± 0.66 mya 

Clades E and F 

8.4% 

2.4 mya ± 0.54 mya 

Between other genera 

Tylomys and Ototylomys 

21.3% 

6.1 mya± 0.90 mya 

Xenomys and Hodomys 

16.1% 

4.6 mya± 0.75 mya 

Baiomys and Scotinomys 

17.2% 

4.9 mya± 0.78 mya 

Reithrodontomys and Isthmomys 

18.0% 

5.1 mya± 0.80 mya 


Discussion 


Phylogenetic analysis of Cytb sequences arranged 
samples of TV. sumichrasti into two major clades (Fig. 
2). The first clade (I) contained samples of TV. sum¬ 
ichrasti from western Mexico and Central America 
together with the sample of Otonyctomys hatti. The 
second clade (II) contained the two samples of TV. 
sumichrasti from southern Mexico. The inclusion 
of the sample of O. hatti with the nine samples of TV. 
sumichrasti produced a paraphyletic arrangement and 
inferred a heretofore unrecognized set of relationships 
among the samples of Nyctomys. In addition, levels of 
genetic divergence, as reflected by the genetic distance 
data (Table 1), document an unusually high level of 
differentiation among clades. Together, these data sug¬ 
gest that a reassessment of the taxonomy of Nyctomys 
relative to Otonyctomys is necessary. Although the 
geneology we recovered represents a gene tree and our 
sample size and geographic coverage is inadequate for 
a formal revision, several comments are appropriate. 


Morphologically, specimens of Nyctomys are 
easily distinguishable from the monotypic Otonycto¬ 
mys. Otonyctomys is smaller in the majority of cranial 
measurements, has much larger auditory bullae and a 
smaller maxillary toothrow (Anthony 1932; Genoways 
et al. 2005); therefore, placing Otonyctomys as a junior 
synonym of Nyctomys is inappropriate. On the other 
hand, the magnitude of genetic differentiation (Table 1) 
between the clades containing samples of Nyctomys and 
Otonyctomys , relative to other closely related genera 
of cricetid rodents, would support recognition of three 
genera: Nyctomys (members of clade II), Otonyctomys , 
and an unnamed genus (members of clade B). For 
example, levels of genetic divergence for comparisons 
of clades I and II (20.0%) and clades A and B (17.7%) 
were comparable for pair-wise values obtained from 
Xenomys/Hodomys , Baiomys/Scotinomys , Reithrodon- 
tomys/Isthmomys , and Tylomys/Ototylomys (ranged 
from 16.1% to 21.3%). Conversion of these distance 
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values into time since divergence estimates would place 
the origin of the unnamed genus in the late Miocene 
(~5.7 ± 0.86 mya), a date comparable to those of other 
Neotomine and Tylomyine genera (Table 1). 

Regardless of generic relationships among clades 
representing Nyctomys and Otonyctomys, multiple spe¬ 
cies likely are included in what currently is recognized 
as N. sumichrasti. Depending on the species concept 
that is invoked (phylogenetic - Cracraft 1983; genetic - 
reviewed by Baker and Bradley 2006), it is possible to 
argue for the recognition of multiple species (perhaps as 
many as five) within the samples included in this study. 
Below, we follow the postulation of Baker and Bradley 
(2006) that genetic distance values >5% (within the 
Cytb gene) exceed the average value for sister species 
of mammals, and as such, should be further evaluated in 
the context that they may represent separate species. 

Clade I is comprised of three minor clades (D, 
E, and F) whose pair-wise distance values range from 
8.4% to 14.4% (Table 1). The sample from Jalisco, 
Mexico (N. s. colimensis , clade D), occurs near sea level 
in tropical deciduous forests, is disjunct geographi¬ 
cally from the remaining samples comprising clades 
E (Guatemala) and F (El Salvador, Honduras, and Ni¬ 
caragua), and differs from them by 14.4% and 13.1%, 
respectively. Samples in clades E ( N . s. decolorus) and 
F (N. s. decolorus and N. s.florencei ) differ by 8.4%. 
This split between low elevation samples of Nyctomys 
is similar to the pattern recovered for Alouatta and Mar¬ 


moset mexicana (Baumgarten and Williamson 2007; 
Gutierrez et al. 2010) and likely involved the Maya 
highlands as a barrier to dispersal between taxa located 
in the Yucatan Peninsula and those found further south. 
Clade II is comprised of samples from Chiapas, Mexico 
(N. s. salvini ) and Oaxaca, Mexico (N. s. sumichrasti ), 
that differ by 10.2%. Both samples are from cloud 
forests, separated by the Isthmus of Tehuantepec. The 
Isthmus has been recognized as an effective vicariant 
barrier for other mid- to high-elevation rodent taxa (Sul¬ 
livan et al. 1997, 2000; Carleton et al. 2002; Edwards 
and Bradley 2002; Arellano et al. 2005; Leon-Paniagua 
2007; Rogers et al. 2007). 

Given the complex geography of southern Mex¬ 
ico and Central America (see Almendra and Rogers 
2012 for a recent summary), it is possible that multiple 
lowland and highland forms exist, especially given the 
influence of the many isolated mountain ranges. Also, 
given the history of the Isthmus of Tehuantepec that 
included climatic changes coupled with marine incur¬ 
sions (Beard et al. 1982; Toledo 1982), it is possible 
that multiple invasions (from north to south, or the 
reciprocal) have taken place. Resolving these issues 
will require thorough geographic sampling and should 
include the addition of samples representing a more 
thorough geographic coverage. In addition, inclusion 
of unrepresented subspecies, ideally material from 
type localities, is essential for determining the prior¬ 
ity of available names in cases where taxa should be 
elevated to species. 


Acknowledgments 


Fieldwork was funded in part by J. Sowell 
through the Natural Science Research Laboratory, 
Museum of Texas Tech University, and the Monte L. 
Bean Life Science Museum and the Office of Research 
and Creative Activities, Brigham Young University. 
Collecting permits were provided by the Consejo Na- 
cional de Areas Protegidas (CONAP), Guatemala, to 
N. Ordonez-Garza. Permits for fieldwork in Mexico 
were issued by the Secretaria de Medio Ambiente, 


Recursos Naturales y Pesca (SEMARNAP) to F. X. 
Gonzalez-Cozatl and F. A. Cervantes. We thank A. 
Almendra-Villalba, E. Arellano, A. Ferguson, F. X. 
Gonzalez-Cozatl, R. Mercado, L. Siles for assisting 
with field work, as well as H. Garner, K. McDonald, 
R. J. Baker, and the staff at the NSRL for tissue loans. 
We also thank C. W. Thompson, E. K. Roberts, H. R. 
Huynh, and M. R. Mauldin for comments on earlier 
versions of this manuscript. 


Corley et al.—Paraphyly in Nyctomys 


7 


Literature Cited 


Almendra, A. L., and D. S. Rogers. 2012. Biogeography 
of Central American mammals: Patterns and pro¬ 
cesses. Pp. 203-229 in Bones, Clones and Biomes: 
The History and Geography of Recent Neotropical 
Mammals (B. D. Patterson and L. P. Costa, eds.). 
University of Chicago Press, Illinois. 

Anthony, H. E. 1932. A new genus of rodents from Yucatan. 
American Museum Novitates 586:1-3. 

Arbogast, B. S., and J. B. Slowinski. 1998. Pleistocene 
speciation and the mitochondrial DNA clock. Sci¬ 
ence 282:1955a. 

Arellano E., F. X. Gonzalez-Cozatl, and D. S. Rogers. 2005. 
Molecular systematics of Middle American harvest 
mice Reithrodontomys (Muridae), estimated from 
mitochondrial Cytochrome b gene sequences. Mo¬ 
lecular Phylogenetics and Evolution 37:52-540. 

Baker, R. J., and R. D. Bradley. 2006. Speciation in mam¬ 
mals and the genetic species concept. Journal of 
Mammalogy 87:643-662. 

Baumgarten, A., and G. B. Williamson. 2007. The distri¬ 
butions of howling monkeys ( Alouatta pigra and 
A. palliata) in southeastern Mexico and Central 
America. Primates 48:310-315. 

Beard, J. H., J. B. Sangree, and L. A. Smith. 1982. Qua¬ 
ternary chronology, paleoclimate, depositional se¬ 
quences, andeustatic cycles. American Association 
of Petroleum Geologists Bulletin 66:158-169. 

Bradley, R. D., and J. Ensink. 1987. Karyotypes of five 
cricetid rodents from Honduras. Texas Journal of 
Science 39:171-175. 

Bradley, R. D., C. W. Edwards, D. S. Carroll, and C. W. 
Kilpatrick. 2004. Phylogenetic relationships of 
Neotomine-Peromyscine rodents: based on DNA 
sequences from the mitochondrial cytochrome b 
gene. Journal of Mammalogy 85:389-395. 

Carleton, M. D., O. Sanchez, and G. Urbano-Vi dales. 2002. 
A new species of Habromys (Muroidea: Noetomi- 
nae) from Mexico, with generic review of species 
definitions and remarks on diversity patterns among 
Mesoamerican small mammals restricted to humid 
montane forests. Proceedings of the Biological 
Society of Washington 115:488-533. 

Cervantes, F. A., J. Nahu Ramirez-Vite, S. Ramirez-Vite, 
and C. Ballesteros. 2004. New records of mam¬ 
mals from Hidalgo and Guerrero, Mexico. The 
Southwestern Naturalist 49:122-124. 


Cracraft, J. 1983. Species concepts and speciation analyses. 
Pp. 159-187 in Current Ornithology (R. F. Johnston, 
ed.). Vol. 1. Plenum Press, New York. 

De Saussure, M. H. 1860. Note sur quelques mammiferes 
du Mexique par M. H. De Saussure. Revue Et 
Magasin de Zoologie Pure Et Appliquee. Series 2 
Volume 12: 97-110. 

Edwards C.W., and R. D. Bradley. 2002. Molecular sys¬ 
tematics of the genus Neotoma. Molecular Phylo¬ 
genetics and Evolution 25:489-500. 

Emmons, L. H. 1997. Neotropical rainforest mammals: A 
field guide. University of Chicago Press, Illinois. 

Genoways, H. H., R. M. Timm, and M. D. Engstrom. 2005. 
Natural history and karyology of the Yucatan 
vesper mouse, Otonyctomys hatti. Pp. 213-218 in 
Contribuciones Mastozoologicas en Homenaje a 
Bernardo Villa-R (V. Sanchez-Cordero and R. A. 
Medellin, eds.). Instituto de Biologia e Instituto de 
Ecologia, UN AM; CONABIO, Mexico. 

Gutierrez, E. E., S. A. Jansa, and R. S. Voss. 2010. Mo¬ 
lecular systematics of mouse opossums (Didel- 
phidae: Marmosa): Assessing species limits using 
mitochondrial DNA sequences, with comments 
on phylogenetic relationships and biogeography. 
American Museum Novitates 3692:1-22. 

Haiduk, M. W., C. Sanchez-Hernandez, and R. J. Baker. 

1988. Phylogenetic relationships of Nyctomys and 
Xenomys to other cricetine genera from G-banded 
chromosomes. The Southwestern Naturalist 
33:397-403. 

Hall, E. R. 1981. The mammals of North America. John 
Wiley and Sons, New York. 

Huelsenbeck, J. P, and F. Ronquist. 2001. MrBayes: 
Bayesian inference of phylogeny. Bioinformatics 
17:754-755. 

Hunt, J. L., J. E. Morris, and T. L. Best. 2004. Nyctomys 
sumichrasti. Mammalian Species 754:1-6. 

Irwin, D. M., T. D. Kocher, and A. C. Wilson. 1991. Evolu¬ 
tion of the cytochrome b gene in mammals. Journal 
of Molecular Evolution 2:37-55. 

Kimura, M. 1980. A simple method for estimating evolution¬ 
ary rate of base substitutions through comparative 
studies of nucleotide sequences. Journal of Mo¬ 
lecular Evolution 16:111-120. 


8 


Occasional Papers, Museum of Texas Tech University 


Lee, M. R., and F. F. B. Elder. 1977. Karyotypes of eight 
species of Mexican rodents (Muridae). Journal of 
Mammalogy 58:479-487. 

Leon-Paniagua, L., G. Navarro-Sigiienza, B. E. Hernandez- 
Banos, and J. C. Morales. 2007. Diversification of 
the arboreal mice of the genus Habromys (Rodentia: 
Cricetidae: Neotominae) in the Mesoamerican 
highlands. Molecular Phylogenetics and Evolu¬ 
tion 42:653-664. 

Miller, J. R., and M. D. Engstrom. 2008. The relation¬ 
ships of major lineages within peromyscine ro¬ 
dents: A molecular phylogenetic hypothesis and 
systematic reappraisal. Journal of Mammalogy 
89:1279-1295. 

Musser, G. M., and M. D. Carleton. 2005. Superfamily 
Muroidea. Pp. 894-1531 in Mammal Species of 
the World: A Taxonomic and Geographic Refer¬ 
ence (D. E. Wilson and D. M. Reeder (eds.). Third 
edition. The Johns Hopkins University Press, 
Baltimore, Maryland. 

Palumbi, S. R. 1996. Nucleic acids I: The polymerase chain 
reaction. Pp. 205-247 in Molecular Systematics (D. 
M. Hillis, C. Moritz, and B. K. Mable, eds.). Sec¬ 
ond edition. Sinauer, Sunderland, Massachusetts. 

Peppers, L. L., and R. D. Bradley. 2000. Cryptic species 
in Sigmodon hispidus : Evidence from DNA se¬ 
quences. Journal of Mammalogy 81:332-343. 

Peppers, L. L., D. S. Carroll, and R. D. Bradley. 2002. 
Molecular systematics of the genus Sigmodon 
(Rodentia: Muridae): Evidence from the mitochon¬ 
drial cytochrome b gene. Journal of Mammalogy 
83:396-407. 

Posada, D., and K. A. Crandall. 1998. MODELTEST: test¬ 
ing the model of DNA substitution. Bioinformatics 
14:817-818. 

Reeder, S. A., and R. D. Bradley. 2004. Molecular systemat¬ 
ics of neotomine-peromyscine rodents based on the 
dentin matrix protein 1 gene. Journal of Mammal¬ 
ogy 85:1194-1200. 

Reeder, S. A., D. S. Carroll, C. W. Edwards, C. W. Kilpatrick, 
and R. D. Bradley. 2006. Neotomine-peromyscine 
rodent systematics based on combined analyses of 
nuclear and mitochondrial DNA sequences. Mo¬ 
lecular Phylogenetics and Evolution 40:251-258. 

Reeder, S. A., and R. D. Bradley. 2007. Phylogenetic rela¬ 
tionships of neotomine-peromyscine rodents using 
DNA sequences from intron 7 of the beta fibrinogen 
gene. Pp. 883-900 in The Quintessential Naturalist: 
Honoring the Life and Legacy of Oliver P. Pearson 
(D. A. Kelt, E. P. Lessa, J. A. Salazar-Bravo, and J. 


L. Patton, eds.). University of California Publica¬ 
tions in Zoology, Berkeley. 

Reid, F. A. 2009. A field guide to the mammals of Central 
America and Southeast Mexico. Second edition. 
Oxford University Press, New York. 

Rogers, D. S., C. C. Funk, J. R. Miller, and M. D. Engstrom. 
2007. Molecular phylogenetic relationships among 
crested-tailed mice (genus Habromys). Journal of 
Mammalian Evolution 14:37-55. 

Saiki, R. K., et al. 1988. Primer-directed enzymatic am¬ 
plification of DNA with a thermostable DNA 
polymerase. Science 239:487-491. 

Sanchez-Hemandez, C., M. de Lourdes Romero-Almaraz, 
R. D. Owen, A. Nunez-Garduno, and R. Lopez- 
Wilchis. 1999. Noteworthy records of mammals 
from Michoacan, Mexico. The Southwestern 
Naturalist 44:231-235. 

Smith, M. F., and J. L. Patton. 1993. The diversifica¬ 
tion of South American rodents: Evidence from 
mitochondrial sequence data for the akodontine 
tribe. Biological Journal of the Linnean Society 
50:149-177. 

Sullivan, J., A. Markert, and C. W. Kilpatrick. 1997. Phy- 
logeography and molecular systematics of the 
Peromyscus aztecus species group (Rodentia: 
Muridae) inferred using parsimony and likelihood. 
Systematic Biology 46:426-440. 

Sullivan, J., E. Arellano, and D. S. Rogers. 2000. Compara¬ 
tive phylogeography of Mesoamerican highland 
rodents: Concerted versus independent response to 
past climatic fluctuations. The American Naturalist 
155:755-768. 

The Mathworks, Inc. 2006. MATLAB® version 6.5. http:// 
www. mathworks. com/ company. 

Tiemann-Boege, I., C. W. Kilpatrick, D. J. Schmidly, and 
R. D. Bradley. 2000. Molecular phylogenetics of 
the Peromyscus boylii species group (Rodentia: 
Muridae) based on mitochondrial cytochrome b 
sequences. Molecular Phylogenetics and Evolu¬ 
tion 16:366-378. 

Timm, R. M., and R. K. LaVal. 2000. Mammals. Pp. 223- 
244 in Monteverde: Ecology and Conservation of 
a Tropical Cloud Forest (N. M. Nadkarni and N. 
T. Wheelwright, eds.). Oxford University Press, 
New York. 

Toledo, V. M. 1982. Pleistocene changes of vegetation in 
tropical Mexico. Pp. 93-111 in Biological Diversi¬ 
fication in the Tropics (G. T. Prance, ed.). Columbia 
University Press, New York. 


Corley et al.—Paraphyly in Nyctomys 


9 


Whiting, A. S., A. M. Bauer, and J. W. Sites, Jr. 2003. Phy- 
logenetic relationships and limb loss in sub-Saharan 


African scincines (Squamata: Scincidae). Molecu¬ 
lar Phylogenetics and Evolution 29:582-598. 


Addresses of authors : 

Megan S. Corley 

Department of Biological Sciences 
Texas Tech University 
Lubbock, Texas 79409-3131 
megan. corley@ttu. edu 

Nicte Ordonez-Garza 

Department of Biological Sciences 
Texas Tech University 
Lubbock, Texas 79409-3131 
nicte. ordonez-garza@ttn. edu 


Duke S. Rogers 

Department of Biology and 
Monte L. Bean Life Science Museum 
Brigham Young University 
Provo, Utah 84602 
duke_rogers@byu. edu 

Robert D. Bradley 

Department of Biological Sciences and 
Natural Science Research Laboratory, The Museum 
Texas Tech University 
Lubbock, Texas 79409-3131 
robert. bradley@ttu. edu 


Editor for this manuscript was Robert J. Baker. 


10 


Occasional Papers, Museum of Texas Tech University 


Appendix 

Specimens from which the cytochrome-Z? gene was sequenced are listed below with: locality informa¬ 
tion, museum voucher numbers, collection or tissue number, and GenBank accession numbers. The remaining 
sequences were obtained from GenBank (sequences generated in previous studies) and are listed by GenBank 
accession number. 

Nyctomys sumichrasti colimensis. —MEXICO: Jalisco; 6 km SE Chamela, UNAM Estacion de Biologia 
(Locality 1, TTU37743, TK19590, JQ183066). 

Nyctomys sumichrasti decolorus .—GUATEMALA: Peten; San Jose; Limite Oeste Biotopo San Miguel 
La Palotada-El Zotz, 17.17241°N, 89.88585°W, elev = 215 m (Locality 3, TTU115406, TK169284, JN851816). 
HONDURAS: Atlantida; Lancetilla Botanical Garden (Locality 4, TTU84484, TK101827, AF195801). 

Nyctomys sumichrasti florencei. —EL SALVAOR: La Paz; about 3 mi NW San Luis Talpa (Locality 5, 
TTU111528,TK34799, JQ183063;TTU111529, TK34801,JQ183064). NICARAGUA: Chinandega; Chichigalpa, 
Belle Vista (Locality 6, TTU105094, TK113517, JQ183061); Atlantico Norte; Siuna, El Balsamo (Locality 7, 
TTU100336, TK121424, JQ183062). 

Nyctomys sumichrasti salvini. —MEXICO: Chiapas; Municipio Chamela, Cerro Tzontehuitz, 13 km NE 
San Cristobal de las Casas, 2,880 m (Locality 8, BYU14466, DSR4226, JQ183065). 

Nyctomys sumichrasti sumichrasti .—MEXICO: Oaxaca; Distrito Ixtlan, 28 km S W (by road) La Esperanza, 
17°35’08”N, 96°30’41”W, 2,950 m (Locality 9, CMC101, DSR5788, JQ183067). 

Reference samples obtained from GenBank: 

Neotominae.— Baiomys taylori (AF548469), Scotinomys teguina (JN851815), Hodomys alleni (AF186801), 
Xenomys nelsoni (AF307838), Neotoma mexicana (AF294345), Ochrotomys nutalli (AY195798), Onychomys 
arenicola (AY195793), Reithrodontomys fulvescens (AF176257), Isthmomys pirrensis (FJ214681), Habromys 
ixtlani (DQ973099), Neotomodon alstoni (AY195796), Megadontomys thomasi (AY195795), Podomysfloridanus 
(EF989977), Osgoodomys banderanus (DQ000473), Peromyscus californicus (AF155393), and Peromyscus 
maniculatus (AY322503). 

Tylomyinae. —Otonyctomys hatti (JQ183060), Tylomys nudicaudus (AF307839), and Ototylomysphyllotis 
(AY009789). 

Sigmodontinae. —Oiyzomyspalustris (DQ185382) and Sigmodon hispidus (AF155420). 


